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Secondary antiproton flux induced by cosmic ray interactions with the atmosphere 
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The atmospheric secondary antiproton flux is studied for detection altitudes extending from sea 
level up to about 3 earth radii, by means of a 3-dimensional Monte-Carlo simulation, successfully 
applied previously on other satellite and balloon data. The calculated antiproton flux at mountain 
altitude is found in fair agreement with the recent BESS measurements. The flux obtained at 
balloon altitude is also in agreement with calculations performed in previous studies and used for 
the analysis of balloon data. The flux at sea level is found to be significant. The antineutron flux 
is also evaluated. The antiproton flux is prospectively explored up to around 2 10 4 km. The results 
are discussed in the context of the forthcoming measurements by large acceptance experiments. 

PACS numbers: 94.30.Hn,95.85.Ry,96.40.-z,13.85.-t 



I. INTRODUCTION 

The antiproton (p) has a particular status in the spec- 
trum of Cosmic Radiation mainly because of its particu- 
lar production dynamics and kinematics. The main part 
of the Cosmic Ray (CR) p spectrum measured in balloon 
and satellite experiments is well accounted for by assum- 
ing it to consist of a secondary flux, originating from the 
interaction between the nuclear CR flux and the inter- 
stellar matter in the galaxy (ISM) pj. It is expected to 
be dominant with respect to components from other pos- 
sible origins. Such other contributions of primary origin 
and of major astrophysical interest, have been considered 
recently. In particular, the p flux induced by annihila- 
tion of dark matter constituents Jp, 0, Q , and by primor- 
dial black hole evaporation 0, [|| have been discussed. 
All these possible contributions are intimately entangled 
together and their phcnomenological discntagling relies 
critically on the accuracy of the experimental data. The 
measurements of the p flux thus provide a sensitive test of 
the production source and mechanism, and of the prop- 
agation conditions in the galaxy 0, 0, S 13 fill E3 ■ 

CR antiprotons have been experimentally studied for 
several decades by satellite or balloon borne experiments 
(see references in Eg). Several recent balloon experi- 
ments, like BESS O E3 and CAPRICE have 
collected new data samples whose analysis have provided 
determinations of the galactic p flux over a kinetic energy 
range extending from about 0.2 GeV kinetic energy up to 
about 50 GeV. In these works, the values of the antipro- 
ton galactic flux were obtained by subtracting the calcu- 
lated atmospheric p flux from the values of the measured 
total flux. 

Secondary galactic as well as atmospheric antiprotons 
are both produced in hadronic collisions by the same el- 



ementary reaction mechanism in nucleon-nucleon colli- 
sions between the incident CR flux and either ISM nu- 
clei (mainly Hydrogen) in the galaxy, or atmospheric 
nuclei (mainly Nitrogen) in the atmosphere. The basic 
p production reaction is the inclusive NN — > pX pro- 
cess, N standing for nucleon and X for any final quan- 
tal hadronic state allowed in the process. The ratio of 
p production in the galaxy and in the atmosphere scales 
with the ratio of matter thickness (in units of interaction 
length Xj) crossed by protons in the two media. These 
thicknesses are known to be of the same order of magni- 
tude. In addition, both flux are driven by similar trans- 
port equations (see |t1 IsL ITtI ITsj for example), with how- 
ever the escape term arising from convective and diffusive 
effects in the interstellar medium on the galactic flux Q , 
making a significant difference with the transport of flux 
in atmosphere, which tends to decrease the transported 
flux compared to the atmospheric transport conditions. 

It can be shown using a Leaky Box Model (LBM) for 
the galactic transport and a simple slab model for the 
production in the upper atmosphere |l9j |. that the ratio 
of the atmospheric p flux at balloon altitude N atm to the 
galactic p flux at TOA N ga i, is approximately: 
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Where XAtm is the thickness of atmosphere on top of 
balloon experiments, while A e is the LBM escape length, 
a(pp — > pX) and a(p Atm — ► pX) being the inclusive 
antiproton production cross sections on hydrogen and on 
the atmospheric nuclei respectively, while m# and mAtm 
are the hydrogen and the mean atmospheric nuclear mass 
respectively. Using x y i tm =3.9 g/cm 2 (for a 38 km alti- 
tude) [2(J and A e =8 g/cm 2 and 11.8 g/cm 2 for particle 
rigidities of 3 GV and 10 GV respectively |2l|, the above 
ratio is found to be of the order of 0.15 and 0.2, respec- 
tively. The contribution of the atmospheric antiproton 
production to the total flux measured in balloon experi- 
ments is thus not expected to be negligible with respect 
to the galactic component. The correction of the total 
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flux from the atmospheric contribution therefore needs 
the latter to be calculated very carefully since the accu- 
racy on the evaluation of this component sets a limit of 
accuracy on the final value of the measured galactic flux. 

It must be emphasized that studying also the sec- 
ondary proton flux in the atmosphere in this context 
is interesting since the latter is very sensitive to all the 
components of the simulation process, in particular to 
the secondary proton production cross section which con- 
tributes as well to the generation of the antiproton flux. 
The comparison of the calculation to the recently mea- 
sured data provides a robust validation of the approach 
used and of the overall calculation, and firmly supports 
the reliability of the results reported here. This study is 
reported in a separate (companion) paper |22^ . 

The present work is an extension of a resarch program 
aiming at the interpretation of satellite data and which 
first results on the flux of protons, leptons, and light ions, 
below the geomagnetic cutoff (GC) at satellite altitude, 
have been reported recently |3 H3| • 

The paper reports on the calculated p atmospheric flux 
over the range from sea level up to satellite altitudes by 
Monte-Carlo simulation. The main features of the calcu- 
lations are described in section [HI The production cross 
sections used in the event generator are given in section 
IIIII The results are discussed in section llVl The work is 
concluded in section IV! 



II. SIMULATION CONDITIONS 

As mentioned above, the flux of secondary atmospheric 
antiprotons has been investigated using the same simula- 
tion approach which has allowed to successfully account 
for the p, d, He, and e ± experimental flux below the Geo- 
magnetic Cutoff (GC) measured by the AMS experiment, 
as well as the experimental proton and muon flux in the 
atmosphere, the latter being studied together with the 
atmospheric neutrino flux |26j . 

The same computing environment has been used here 
for the charged particle propagation in the terrestrial en- 
vironment including the atmosphere, as in the previous 
studies, with the event generator being dedicated how- 
ever, based on the antiproton production cross section in 
nucleon-nucleon collisions. 

The calculation proceeds by means of a full 3D- 
simulation program. Incident Cosmic Rays are gener- 
ated on a virtual sphere chosen at a 2000 km altitude. 
Random events are generated uniformly on this sphere. 
The local zenith angle distribution of the particule mo- 
mentum is proportional to cas(9 z )dcos(6 x ), 9 Z being the 
zenithal angle of the particle, in order to get an isotropic 
flux at any point inside the volume of the virtual sphere. 
The geomagnetic cut-off is applied by back-tracing the 
particle trajectory in the geomagnetic field, and keeping 
in the sample only those particles reaching a backtracing 
distance of 10 Earth radii. Flux conservation along any 
allowed particle path in the geomagnetic field is ensured 



by application of Liouville's theorem. The normal parti- 
cle propagation as well as its back-tracing are performed 
using the adaptative R ung e-Kutta integration method in 
the Geomagnetic field [2a ]. 

1. For the incident CR proton and helium flux, func- 
tional forms fitted to the 1998 AMS measurements were 
used [H H| (see also [H IM Elj). The heavier com- 
ponents of the CR flux were not taken into account in 
the calculations (see |2(j). For other periods of the solar 
cycle, the incident cosmic flux are corrected for the dif- 
ferent solar modulation effects using a simple force law 
approximation 32] . 

2. Each particle is propagated in the geomagnetic field 
and interacts with nuclei of the local atmospheric den- 
sity according to their total reaction cross section and 
producing secondary nucleons p, n, and antinucleons p , 
n , with cross sections and multiplicities. This important 
issue is discussed in section ITTT1 below. The specific ion- 
ization energy loss is computed for each step along the 
trajectory. 

3. In the following step, each secondary particle pro- 
duced in a primary collision is propagated in the same 
conditions as incident CRs in the previous step, resulting 
in a more or less extended cascade of collisions through 
the atmosphere, which may include up to ten generations 
of secondaries for protons for the sample generated in this 

work mm. 

For the antinucleon inelastic collisions, only the anni- 
hilation reaction channel was taken into account. Non 
annihilating inelastic N + A — > N + X (N standing for 
antinucleon) interactions whose contribution to the total 
reaction cross section <tr, is small. It consists basically 
of the single diffractive dissociation cross section (for the 
proton target in individual Np collisions), and it would 
be of the same order of magnitude as for pp colisions, 
namely ~10% of an or less at the energies considered 
here [33]. It has been neglected at this stage. It will be 
included in the further developments of the calculation 
program. 

The reaction products are counted whenever they 
cross, upwards or downwards, the virtual detection 
spheres (several can be defined in the program) at the al- 
titude of the detectors: from sea level up to about 36 km 
for ground and balloon experiments (BESS, CAPRICE), 
370 km for the AMS satellite experiment. Higher alti- 
tudes up to more than 10000 km were also investigated, 
with the purpose of understanding the dynamics of the 
population of quasi-trapped particles in the earth envi- 
ronment (see section flV C|) . Each particle is propagated 
until it disappears by nuclear collision (annihilation for 
antinucleons), stopping in the atmosphere by energy loss, 
or escaping to the outer space beyond twice the genera- 
tion altitude I2II I2H Eif . 

In the terminology used in the following, one event is 
defined as the full cascade induced by an incoming CR 
particle interacting with one atmospheric nucleus. For 
each CR producing at least one secondary particle, the 
whole event is stored with all the relevant topological, 
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dynamical, kinematical, and geographical informations. 
This includes the collision rank, geophysical location, al- 
titude, momentum and particle type, and parent particle 
type, in form of event files. The collision rank is defined 
as the number of a given collision in the cascade initiated 
by the first CR interaction with atmosphere (rank 1). 

The calculations do not include any adjustable param- 
eter. 



III. CROSS SECTIONS 
A. Proton induced secondaries 

1. Protons 

The inclusive p + A — ► p + X proton production cross 
sections used are described in ref |23j . They are based on 
the results of refs [34| an d HH f° r the two components 
corresponding basically to forward (or direct quasi elas- 
tic) and backward (or relaxed deep inelastic) productions 
respectively. The values obtained have been found in rea- 
sonable agreement with the results of the INCL model of 
intranuclear cascade calculations |3fij . This cross section 
allows to reproduce very successfully the atmospheric sec- 
ondary proton flux down to the lowest altitude p3 |. 

The cross section for secondary neutron production 
was taken equal to that of proton production. Similarly, 
the neutron induced cross sections of secondary nucleon 
production used was taken the same as for protons, the 
coulomb interaction making negligible differences over 
the considered energy range. 



TABLE I: Values of the parameters of relation^ obtained in 
fitting the data of figure [71 



Parameter Ci Ci Cz Ca C5 Ce \x bo a 
value 0.042 5.92 0.96 2.19 84.3 10.5 1.1 0.12 2.24 



With b{pt) = b Q ■ p t within the range of 4- momentum 
transfer considered here. Figure Q shows some of the 
results obtained ^3 by fitting the set of available data 
studied for this work with relation The energy de- 
pendent exponential factor in the second term of relation 
El was introduced for it was found that this term con- 
tributes only at low incident energy. The values of the 
parameters obtained are given in table [I] These values 
are si gnifi cantly different from those given in the original 
work |42t l43j on a smaller number of data. This work is 
currently being extended up to 400 GeV incident proton 
energy [44| . 



B. 4 He induced secondaries 

1. Protons 

The inclusive A He + A — > p(n) + X proton (neutron) 
production cross section was obtained as described in sec- 
tion IIII A II above for the p + A — > p + X cross section 
(using the total reaction cross section from for this 
system), renormalized to the available ex per imental mul- 
tiplicities measured for this reaction jig, [43 • 



2. Antiprotons 



The inclusive p + A — > p + X antiproton production 
cross section has been obtained by fitting a set of avail- 
able experimental data between 12 GeV incident kinetic 
energy and 24 GeV/c incident momentum |3^ . l38ll3^ . l40| . 
using a modified version of the analytical formula pro- 
posed in ref |4l| 1 the latter being referred to as KMN 
(for Kalinovski, Mokhov, Nikitin) in the following. 

The invariant triple differential cross section is de- 
scribed by means of the formula used in |4l| : 



(F d3<T ) 



O-R 



Cii 6fe '(l - x) c - exp(-C 3 x)0>(p t ) (1) 



In this relation the kinematical variables pt and x are 
the transverse momentum and the fractional energy of 
the particle respectively as defined in [41j (relation 3.26 
in this reference), while ctr is the total p + A reaction 
cross section. The function $(pt) was modified as: 



$(pt) = exp(-C 4 p?; 



exp(-C 6 x t ) r 
5 (p^ + /i 2 ) 4 exp( ~ Q%/s) (2) 



2. Antiprotons 

The inclusive antiproton production cross section was 
evaluated by means of the wounded nucleon model 
[48L I49I ] using the experimental values of the total re- 
action cross sections for the 4 He + A and p + A systems, 
and of the p production multiplicity in nucleon-nucleon 
collisions [iij available. 

In this model, the particle production multiplicity 

< tiab > in the collision between ions A and B is re- 
lated to the multiplicity nucleon-nucleon (AfAf) collision 

< nj\fj^ > by the relation: 



< n A B > 



A 



0~ P B 
0~AB 



B 



Q~pA 
0~AB 



< > (3) 



with cry being the total reaction cross section between 
the i and j system. Using this model, the p production 
multiplicity induced by the CR He component on the 
nitrogen component of the atmosphere is found to be 
< n-HeN >~ 1-55 < n P N >~ 2.5 < npp >. 
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FIG. 1: Sample of the fits to the p + A — *• p + X cross section data using the parametrized relation^ given in the text, with 
the parameters of table 171 Top left: Differential cross sections measured at 5.1° forp + C (full squares) and for p + Al 

(full circles) as a function of the particle momentum, compared with fit using^\ (thick solid for Al and thick dashed for C ), 
and using the KMN relation and parameters from fjjj (thin solid and thin dashed resp.). Top right: Same comparison for 14-6 
GeV/c p + Al invariant differential cross section versus mass transfer from \3S) : Fit with^\ (thick line) and KMN calculations 
(thin line). Bottom left: Same for 19.2 GeV/c p + p (full squares) andp + Al (full circles) rapidity distributions from \39\1: 
fit with^\ (solid and dashed thick lines resp.) and KMN calculations TTil (thin solid and dased lines resp.). Bottom right: 24 
Ge V/ c p + Al invariant differential cross sections at various angles ( in msr on the figure ) from ]A&j compared to the fits withU\ 
For each measurement angle above the first (17 mrad), each next cross section has been multiplied by i0 _1 > _2 >-- for presentation 
purpose. KMN calculations are not shown on this figure for legibility. Note that the same definitions of cross sections have 
been used as in the original references. There is a clear target mass dependence of the differential cross section in the top left 
panel, while in the bottom left panel almost no such dependence is observed, because the observable displayed is a multiplicity, 
i.e., ratio of differential cross section to total reaction cross section. 



C. Total reaction cross sections 

• Protons: The values of the proton total reac- 
tion cross section on nuclei used were obtained from the 
parametrization of |50| . and checked on the carbon data 
from |45| . 

• 4 He: The 4 He+A total reaction cross sections used 
were taken from 1451. 



• Antinucleons: The p total reaction cross section was 
taken from |5l| . with the energy dependence from the 
data compilation of [S^ • The same production cross sec- 
tion and total reaction cross section have been assumed 
for n production as for p . 
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IV. RESULTS 

A sample of about 35 10 6 CRs have been simulated, of 
which 20% were effectively propagated to the atmosphere 
(above GC), for detection altitudes going from sea level 
up to 10 4 km altitude, including the BESS/CAPRICE 
balloon altitude (ss38 km), the AMS orbit altitude 
(«370 km), and the recent BESS measurement terrestrial 
altitude (2770 m) . The flux at sea level was calculated to 
investigate the possibility of ground level measurement of 
atmospheric antiprotons with existing devices [53j. This 
was achieved independently by BESS at mountain alti- 
tude and the results are discussed below. 



A. Particle trajectories in the Earth magnetic field 

The time of confinement of particles in the earth en- 
vironment together with their particular trajectories, 
determine their status with respect to the three cate- 
gories: trapped, semi-trapped, and non trapped (escape). 
Trapped particles are spiralling back and forth around 
and along the magnetic field lines long enough to drift 
many times around the earth (see for example [H3. l55j 
and below). Trapped particles are practically not ob- 
served in the energy domain considered here. They are 
not dynamically forbiden however and a few trajectories 
with a few 10 2 bounces are observed, which corresponds 
to short-lived trapped particles. Quasi trapped particles 
are in similar kinematic conditions but accomplish only 
a limited number of bounces at mirror points before be- 
ing absorbed or before escaping (see examples below). 
This concept appeared during the first years of radiation 
belts studies [Hg (see also the discussion in JMD- Es- 
cape particles do not match the kinematic conditions for 
being trapped at their production point and escape in 
a very short time to the outer space. All intermediate 
situations between the stereotypes of quasi-trapped and 
escape trajectories are in fact observed in the simulation 
results (see example in bottom left fig|21)- 

Figure [21 shows four examples of characteristic trajec- 
tories of antiprotons generated in this study. Each of the 
4 panels gives a side view (projection on the meridian 
plane, top left), side view zoomed around the production 
point showing the spiralling trajectory of the particle (top 
right), top view (projection on the equatorial plane, bot- 
tom left), and 3-D representation (bottom right) of each 
of the selected trajectories. 

The top left event (1.52 GeV kinetic energy) is an 
escape particle produced close to the North pole. Top 
right is a semi trapped single bounce event (0.54 GeV) 
annihilating in the atmosphere close to its production 
point. Bottom left is a longer lifetime, multi bounces, 
semi trapped event (2.48 GeV), drifting around the earth 
about three quarters of a turn before annihilation in the 
atmosphere in the SAA region. Bottom right is inter- 
mediate between semi-trapped (since it displays at least 
one clear bounce) and escape event (0.54 GeV). It is a 



type of event for which the first adiabatic invariant (mag- 
netic moment conservation) is not conserved because of 
a large variation of the magnetic field along the radius of 
gyration [54| . 

B. General features of the simulated data 

Figure [3] shows a few basic distributions of physics 
observables relevant to the dynamics of the process 
for two detection altitudes: 38 km (solid line) and 
380 km (dashed line) corresponding to balloon and satel- 
lite (AMS) altitudes, inside and outside the atmosphere, 
respectively. The rank of the collision producing the an- 
tiproton (top left panel on the figure) appears to extend 
from 1 up to about 10 for the simulated sample, showing 
that ps are produced up to the tenth generation in the 
collisions cascade. The distributions are a little different 
for the two altitudes, with a significantly larger number 
of p occuring from first interaction at the lower altitude. 
The altitude distribution of the production point for the 
detection at 38 km (top right) shows a discontinuity at 
this altitude due to the incoming flux dominated by pro- 
duction from the upper layer of atmosphere. The mean 
production altitude is found around 46 km and 48 km for 
the lower and upper detection altitudes respectively. The 
particle momentum spectrum at the production point 
(bottom left) is found harder at the higher altitude. The 
number of bounces at the mirror points for particle tra- 
jectories spiralling around the magnetic field lines are 
found as expected very different for the two altitudes of 
detection: At 38 km, only a small population is seen to 
reach a number of bounces larger than a few units (5-6) . 
This flux is significant however and must correspond to 
trajectories lying mostly outside of the atmosphere. At 
380 km, the observed flux of p trajectories with more 
than one bounce is larger by about 2 orders of magni- 
tude than at 38 km, corresponding to the population of 
quasi-trapped particles as discussed previously in p3ll58| . 



C. Antiproton flux at mountain altitude 

The recent measurements of the p flux at 2770 m of 
altitude by the BESS collaboration 59] allows a sensitive 
test of the ability of the present simulation program to 
account for the observed flux since at this altitude the 
p production occurs on the average after a casade of 4 
collisions on the average (see previous section and figEJ). 
Another highly sensitive test of the overall calculation is 
provided in [22J on the atmospheric proton flux. 

Figure 01 shows the p spectrum at 2770 m of altitude 
measured by BESS, compared to the simulation results. 
The latter has been run with the geometrical acceptance 
function of the BESS spectrometer given in [5^ (figures 
4.37 to 4.39, see also |6lj), the overall acceptance angle 
being of the order of 25° . The total CR p + He and par- 
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FIG. 2: Examples of antiproton trajectories in the 



tial He flux are shown on the same figure. The He flux 
contribution is seen to produce a small fraction of about 
5% of the full p flux. Although the total flux calculated 
somewhat underestimates the experimental values, the 
overall agreement is good, the calculated values being on 
the average within one standard deviation from the ex- 
perimental values. This gives confidence in the results of 
the calculations obtained for the other altitudes investi- 
gated and reported below. 

Figure[S]compares the experimental zenith angle distri- 
butions of the p flux to the calculated values (histograms) 
for the same kinetic energy bins as measured by BESS 




magnetic field. Details are discussed in the text. 



[59j. On this figure, the overall agreement between data 
and calculations again appears to be good for all energy 
bins. 

Note that no upward particle have been produced at 
this altitude in the simulated sample, as it could be ex- 
pected 0. 



D. Balloon data 

In this section the atmospheric p flux at ballon alti- 
tude is investigated for comparison with the atmospheric 
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FIG. 3: General features of the simulated p sample at ballon 
(38 km) and satellite (380 km) altitude. Top left: Rank dis- 
tribution (see text); Top right: Altitudes of production; Bot- 
tom left: Momentum distributions; Bottom right: Numbers of 
bounces effected by the particle between the mirror points. The 
spikes observed for high bounce multiplicity in this distribu- 
tion correspond to quasi trapped particles trajectories crossing 
many times the detection altitude. 



p corrections made to the raw flux data in the BESS and 
CAPRICE experiments. 

Figure [fj] shows the values of the galactic p flux ob- 
tained from the BESS and CAPRICE measurements re- 
spectively. These values have been obtained from the 
measured raw flux by subtraction of the atmospheric 
p flux evaluated using an avera ge o f theoretical calcu- 
lations for the BESS experiment 15], and using the cal- 
culations of ref H3 for CAPRICE. On the figure, the 
atmospheric flux calculated in |60j is compared with the 
results from the present work (see also 17]). For the 
two sets of data, it appears that the the present calcula- 
tions are in fairly good agreement with the atmospheric 
antiproton flux obtained from transport equation calcu- 
lations and used to correct the raw flux measured. There 
is a slight trend however for the simulation results to be 
larger than those obtained from the differential equation 
approach by about 20% over the range 10-30 GeV. At 
low p energies the opposite trend is observed and the 
simulation results are found significantly below the val- 
ues obtained from the differential equation. One might 
say that the simulation results should be taken with care 




E. (GeV/n) 



FIG. 4: Antiproton flux data at 2770 m measured by BESS 
\5$j (symbols) compared with simulation results (histograms). 
Thick histogram: full calculation; Thin histogram: CR 4 He 
contribution; dotted line: transport equation calculation from 
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FIG. 5: Experimental zenith angle distribution of the antipro- 
ton flux at 2770 m in bins of kinetic energy compared with 
simulation results. 
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FIG. 6: Left : BESS results: Galactic antiproton distributions deduced from the data (full squares), atmospheric antiproton flux 
from ref \6(\] (curve), and corrections applied by the authors of \lS\j to correct the raw flux values for the atmospheric contribution 
in the original work (open crosses), compared with the atmospheric flux obtained in the present work: all (CR proton + 4 He) 
contributions (thick histogram); CR i He contribution (thin histogram). Right : Same for the CAPRICE experiment data fj jj/ 
(data: open squares; corrections applied in fiffi: star symbols). In both panels Qcgm stands for the geomagnetic latitude of the 
measurements. 




FIG. 7: Atmospheric antiproton spectra expected from the present work for the AMS experiment on the ISS at 380 km of altitude 
for 4 bins of latitudes, compared with the AMS, BESS and CAPRICE data JflA.\lA. wMl in the polar region (bottom right). Full 
line: Downward flux; Dashed line: Upward flux. The increasing flux with decreasing latitude is due to the larger ( quasi- )trapped 
particle population closer to the equator (see text). 
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below 1 GeV because of the lack of experimental cross 
sections for low energy p production, and thus of the 
large corresponding uncertainties on the results of the 
simulation over this range, however, it must be noted 
that the calculated cross sections for low p momentum 
should be in principle reliable for the following reason. 
The p distribution is naturally symmetric in the rapid- 
ity space. The fitting function has the same property 
since it depends only on variables matching this symme- 
try. Therefore a good fit to a set of experimental cross 
sections for p rapidities above the center of mass rapid- 
ity Y cm , automatically ensures the right behaviour of the 
calculated values for rapidities below Y cm , i.e., for par- 
ticle momenta in the laboratory frame, because of the 
symmetry law. 

To conclude this section, the atmospheric p flux cal- 
culated in this work confirm the corrections of the raw 
flux values measured in the BESS and CAPRICE ex- 
periments. This result updates and corrects a previous 
preliminary conclusion on the issue |53| recently quoted 
in |o3| ]. 

The contributions of A>4 CR components were 
not included in the calculations, neither were those 
from the non-annihilating inelastic contributions in the 
p propagation through the atmosphere. These contri- 
butions are small however and not likely to change the 
results by more than a few percents 26]. 



E. AMS altitude 

Future satellite experiments in preparation, plan to 
measure the p flux. A reliable knowledge of the atmo- 
spheric p flux at satellite altitudes is therefore highly de- 
sirable for these experiments to step on explored grounds. 
The p flux calculated for the altitude of the AMS orbit 
are presented in this section. 

Figure shows the expected downwards (Secondaries 
and reentrant Albedo, dashed histogram) and upwards 
(Splash Albedo particles, solid histogram) flux of atmo- 
spheric antiprotons at the altitude of AMS for two re- 
gions of lower geomagnetic latitude: equatorial (0 < 
\8 M \ < 0.2 rad), intermediate (0.2 < |0 M | < 0.4 rad), 
and subpolar (0.8 < \9m\ < 0.9 rad). As expected, the 
flux is predicted larger for the lower latitudes than it 
is around the poles, because of the existence of quasi- 
trapped p components at the low and intermediate lat- 
itudes. Note that the simulated flux is surprisingly 
predicted larger downwards than upwards (bottom left 
panel) . This is in fact an effect of the spectrometer ac- 
ceptance (taken to be 30° with respect to zenith), the 
mean angle for upward particle trajectories being 2 ra- 
dians. The overall upward flux is larger than downward 
by a factor of about 2.5 0, Ell This shows that the 
future satellite measurements of antiproton flux at low 
latitudes will have to be corrected from the atmospheric 
contributions and will probably suffer more uncertainties 
than previously thought. 



The lower right panel compares the p data at TOA 
reported by AMS [Hi, BESS [3, and CAPRICE 0, 
to the flux calculated in the polar region where the AMS 
data have been measured. The calculated (downward) 
atmospheric p component is at the percent level of the 
measured flux in the low energy p range, and can be 
considered as negligible at all energies. 

F. p and n flux at terrestrial altitudes 

The p and n flux have been calculated also at sea level 
in order to provide a realistic order of magnitude of these 
flux for general purpose and for ground testing of em- 
barked experiments. 

• Antiprotons : 

The flux of atmospheric antiprotons at sea level has 
been calculated with the same simulation program. Fig- 
ure 0D shows the distributions obtained at sea level (left) 
and at 4000 m. The energy integrated flux is of the or- 
der of 0.4 10~ 3 p s~ 1 m~ 2 sr~ 1 at all latitudes (see fig 
Ej] below). At 4000 m (right panel on the figure), the 
flux raises to about 7 I0~ 3 p s~ 1 m~ 2 sr~ 1 . These values 
are small but large enough for this flux to be measured 
by currently existing large acceptance detectors (BESS, 
CAPRICE), or in a near future by new detectors under 
construction like AMS and PAMELA. 

• Antineutrons : 

Atmospheric secondary antineutrons may also be of 
interest in ground or balloon measurements [64| • Figure|S| 
shows the kinetic energy spectrum of the expected n flux 
at sea level and at 4000 m (this latter altitude being that 
the Cerro-La-Negra observatory in Mexico where some 
experimental measurements of the antineutron flux are 
being considered 

G. Flux dependence on the altitude 

The p flux has been calculated up to altitudes of 2 
10 4 km with the aim of investigating the general features 
of the dynamics and kinematics of the particles in the 
more remote Earth environment than considered in the 
previous sections. 

Figure El shows the altitude dependence of the energy 
integrated upward and downward p flux in bins of lati- 
tude, assuming a geometrical acceptance of 30° for the 
detector. The calculated distributions display two main 
features: 

1) In the atmospheric range of altitudes, a large peak of 
incoming flux centered around 20 km and corresponding 
to atmospheric secondaries, dominates the distribution 
(< 50km, i.e., <TOA), with basically no associated out- 
going flux (inside the quoted acceptance angle). 

2) Above the atmosphere, surprisingly, the calculated up- 
ward and downward flux are found close to each other up 
to fairly high altitudes, namely sslO km, for the low and 
intermediate latitudes (\0m\ ^$ 0.7 rad). This shows that 
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FIG. 8: Simulation results for the anti-protons (thick line) and antineutrons (thin line) flux at sea level (left) and at 4000 meters 
of altitude (right). 
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FIG. 9: Integrated antiproton flux versus detection altitude 
in bins of latitudes between equator and poles, for downward 
(solid line) and upward (dashed line) p flux. 



a population of quasi-trapped, particles should be ob- 
served in this region of space, i.e., up to around five to 
ten thousands kilometers. This is confirmed by the life- 
time of the particle between their production and their 
absorption, and by the number of bounces of the parti- 
cles between the mirror points of their trajectories, which 
extend up to 100 seconds and several 10 2 bounces (see 
Fig. I3J , respectively for the simulation sample produced. 
At higher latitudes with (\9m\ ^ 0.7 rad), the incoming 
flux progressively disappears, and the outgoing flux then 
corresponds to escape particles. 

From these calculations it can be concluded that there 
should exist a significant flux of quasi-trapped particles 
extending approximately over a decade of altitudes, from 
about 50 km (TOA) up to «10 4 km, depending on the 
particle energy and latitude. This flux has been observed 
already by the AMS experiment in the lower part of the 
altitude range (380 km) [23| ■ Note that the issue was dis- 
cussed long ago in a pioneering paper about the electron 
flux [53. See also (23 (companion paper) for a comple- 
mentary discussion. The energy spectrum of this flux 
extends up to around 10 GeV, which is about the up- 
per momentum limit (8.5 GeV) for which particles can 
match the simple geometrical condition that the gyration 
radius is smaller than the mean trajectory radius to the 
upper atmosphere (for equatorial latitude trajectories, at 
the limit of large, close to 7r/2, pitch angles). 

Future embarked experiments should take these fea- 
tures into account, even though, in principle, the accurate 
knowledge of the kinematical conditions of the particle at 
the detection point allows to know whether it is or not 
of atmospheric origin. 
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V. SUMMARY AND CONCLUSION 

In summary, the secondary antiproton flux produced 
by the Cosmic Ray proton and Helium flux on the at- 
mosphere has been calculated by Monte-Carlo simula- 
tion. The flux calculated for the altitude of 2770 m is in 
fair agreement with the recent BESS measurements. At 
sea level, it is small but measurable and could provide a 
natural facility for testing the identification capability of 
existing experiments or of future devices. For balloon al- 
titudes, the calculated flux has been found in agreement 
with the values calculated in previous works. At satellite 
altitudes (380 km) it appears to be negligible compared 



to the CR flux for polar latitudes, and of the same order 
of magnitude as for the high balloon altitudes for equa- 
torial and intermediate latitudes (below the geomagnetic 
cutoff), indicating that it will have to be taken into ac- 
count in future measurements of the galactic antiproton 
flux at similar altitudes. 
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